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ABSTRACT 



Integral field, or 3D, spectroscopy is the technique of 
obtaining spectral information over a two-dimensional, 
hopefully contiguous, field of view. While there is some 
form of astronomical 3D spectroscopy at all wavelengths, 
there has been a rapid increase in interest in optical and 
near-infrared 3D spectroscopy. This has resulted in the 
deployment of a large variety of integral-field spectro- 
graphs on most of the large optical/infrared telescopes. 
The amount of IFU data available in observatory archives 
is large and growing rapidly. The complications of treat- 
ing IFU data as both imaging and spectroscopy make it 
a special challenge for the virtual observatory. This arti- 
cle describes the various techniques of optical and near- 
infrared spectroscopy and some of the general needs and 
issues related to the handling of 3D data by the virtual 
observatory. 

Key words: Integral-field spectroscopy; Instrumentation; 
Virtual Observatory. 



1. INTRODUCTION 



A basic goal of observational astronomy is to obtain the 
most complete, unbiased description of the objects un- 
der study. A complete description of an astronomical ob- 
servation is inherently multi-dimensional, at each posi- 
tion on the sky one can measure intensity, wavelength or 
energy, and up to four polarization (Stokes) parameters. 
Each of these quantities can also vary with time. Since 
almost all astronomical objects have properties that are 
complex — structures on all scales, asymmetries, gra- 
dients, and variability — there is a strong need to ob- 
tain full 2D (spatial) information of all observables. The 
maximizes the information content of the observations 
and provides the most constraints possible for compari- 
son with models and theories. 

However, the nature of light and the technological limi- 
tations of instruments and detectors in general prevent us 
from making the "perfect" observation. Radio and sub- 
mm interferometric observations come the closest to the 



ideal. A map of the UV plane provides intensity, fre- 
quency (energy), and polarization information. The data 
cube does contain time information as well, though usu- 
ally repeated observations are often needed to study vari- 
ability. Many modern 2D X-ray detectors are also inher- 
ently multi-dimensional since the time and energy of each 
photon are recorded. However, obtaining high energy res- 
olution usually means the loss of one spatial dimension. 

Likewise, in the optical and near-infrared (NIR) the avail- 
able technology limits the number of dimensions that can 
be observed simultaneously. Sensitive photon-counting 
systems such as photo-multiplier tubes, avalanche photo- 
diodes, and superconducting tunnel -junction devices (jde] 
Bruijne et al. 2002| are ideal for studying intensity vari- 
ability on short timescales, but they have low energy res- 
olution and are either single-channel devices or small 
arrays and so are inefficient for studying large areas of 
sky or for doing high-resolution spectroscopy. All large- 
format detectors such as photographic plates and CCDs 
are two dimensional and while they are good at counting 
the number of photons they have poor intrinsic energy 
resolution. Therefore, optics must be used to separate dif- 
ferent wavelengths. Simply placing a dispersive element 
in the light path is one valid solution. While this "slitless" 
spectroscopy can be a useful way of obtaining spectral 
information over a full 2D field, it has the disadvantages 
of high background noise in ground-based observations, 
spectral resolution dependent on image quality, and con- 
fusion in crowded regions. Slitless spectroscopy has been 
used extensively to search for emission lines in limited 
wavelength regions (e.g. |Salzer et al. 2000; Douglas et| 



al. 2002 j ). The background and spectral resolution prob- 
lems are often solved by isolating the region of interest 
with a slit; spectra are obtained only along one physical 
dimension and the second dimension of the detector con- 
tains the wavelength information. This is very appropri- 
ate for compact isolated objects, but mapping large areas 
of extended, complicated sources with a slit is time con- 
suming and it can be difficult to register spatially dithered 
spectra. Also, obtaining spatially contiguous spectra is 
usually not possible. 

However, in the last 30 years clever optical/NIR instru- 
ment designers have developed 3D (x,y,X) techniques 
that produce spectra of moderate to high resolutions (R > 
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100) over a contiguous field. Other advantages of these 
techniques are easier acquisitions (no need to center pre- 
cisely in a slit), radial velocities that are not affected by 
slit-centering errors, and no light losses that are other- 
wise unavoidable with narrow slits. These techniques 
were considered rather specialized at first but in the last 
10 years the advantages of the techniques and the suc- 
cesses of the early 3D, or integral-field, spectrographs 
have made them popular and very common. Most current 
and future large optical/NIR telescopes have or will have 
integral-field options in their standard instrument suites. 

Therefore, IFU (integral-field unit) data is becoming 
more and more common in telescope archives Also, large 
surveys such as SAURON and Atlas3D ( |de Zeeuw et al.| 



2002 Krajnovic 2007) are being carried out now and 



even more ambitious projects are being planned (see be- 
low). Thus, it is important for the Virtual Observatory 
(VO) to support integral-field data. This is a challenge 
because it requires a combination of imaging and spec- 
troscopic techniques. However, since IFU data is still rel- 
atively new and many common reduction packages do not 
include good tools for manipulating IFU data, there is an 
opportunity for the VO to set standards and create useful 
tools that will draw people to the VO. 

The remainder of this proceedings will review the dif- 
ferent types of optical/IR spectrographs and the different 
formats of data that are being produced. Then some basic 
requirements for handling and visualizing IFU data in the 
VO will be discussed. 



2. OPTICAL/IR 3D TECHNIQUES 



The two main types of optical/IR 3D instruments are 
spectrometers and spectrographs. Within each type there 
are a variety of related designs. This section will review 
the basic features of the most common variations. Some 
other good recent technical reviews of different 3D tech- 
niques are Bennett ( 2000| l and AlUngton-Smith ( 2006| l. 



2.1. Spectrometers 

The fundamental principle of a spectrometer is that the 
wavelength region of interest is scanned in discreet steps 
and a narrow-band 2D image is taken at each wavelength, 
thus building up a stack of 2D images, or cube. Interfer- 
ence techniques are used to make the adjustable narrow- 
band filters. 

One of the most common types of 3D spectrometer is the 
imaging Fabry-Perot, sometimes called a tunable filter, in 
which the interference between two closely spaced and 
adjustable partially-transmissive plates creates a narrow 
bandpass. There are also a variety of related techniques 
( Bland-Hawthorn 2000 1. The main advantage of the tun- 
able filter is that it can have a very wide field, up to sev- 
eral arcminutes in diameter. However, the spectral reso- 
lution is relatively low, between 100 and 1000, and often 
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Figure 1. A 3D spectrometer result from tlie imaging FTS 
BEAR on the CFT \Cox et al.\ 0O2\ . This shows the 
velocity-integrated H2 emission from the planetary neb- 
ula NGC 7027. Outflow axes identified from the data 
cube are labeled. 



R and wavelength coverage must be traded off against 
each other Examples of this type of instrument are the 
Taurus tunable filter that has been used on the AAT and 
WHT and GRIF behind the adaptive optics system at the 
CFHT ( IClenetet al.j[2002] l. Also, several new tunable fil- 
ters are being commissioned or developed, including the 
Maryland-Magellan Tunable filter (Dressier et al. 2006) , 
FTF2 for the Flamingos 2 NIR MOS spectrograph for 
Gemini South ( |Scott et al.| [2006|, and the OSIRIS in- 
strument on the GTC ( jCastaiieda et al.||2002[ i. 

One other 3D spectrometer is the Imaging Fourier Trans- 
form Spectrometer (FTS) in which a FTS a Michel- 
son interferometer defines the bandpass. These instru- 
ments have smaller fields-of-view, typically less than 1 
arcminute, but the spectral resolution can be as high as 
about 30000. The FTS is sensitive to shot noise from the 
sky so ground-based observations will often use a band- 
pass filter to limit the spectral range. Figure [T] shows an 
example image from the Bear FTS spectrometer on the 
CFHT ( iCoxetal.||2002] i 

The data format of a 3D spectrometer is a regularly grid- 
ded 3D {x, y, A) datacube. Each (x, y) plane is a typical 
image which can be manipulated by most image process- 
ing software. This data is similar to a radio datacube and 
is the easiest to process and visualize with current soft- 
ware. 

The main drawbacks of spectrometers are that the spec- 
tral range is relatively small and, more importantly, that 
the spectra are not created simultaneously. Therefore, 
any changes in image quality, sky transparency, and back- 
ground during the scan will make the spectra difficult to 
process and interpret. 
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2.2. Spectrographs 



Integral field spectrographs (IFSs) are a more spectra- 
based approach to 3D spectroscopy than the spectrom- 
eters. There are three basic flavors of IFS concepts — 
lenslet only, fiber based, and image slicer — each with 
particular advantages and applications (Figure |2]). In all 
the flavors optics are used to sample the focal plane into 
"spaxels" (spatial pixels) and then the light from each is 
dispersed by a standard grating or grism spectrograph. 
Therefore, the spectral resolution is dependent on the rul- 
ing of the grating/grism and the effective size of the en- 
trance aperture. Most IFSs operate with spectral resolu- 
tions between 500 and 5000. An advantage of the IFS is 
that all of the spectra are obtained simultaneously. How- 
ever, since detector area is needed to store all the spectra 
the field-of-view (FOV) has to be sacrificed. The largest 
fields are about an arcminute in diameter. Techniques 
of modifying the field size include using fore-optics to 
change the scale (arcsec/mm), trading off FOV for spec- 
tral coverage, and by increasing the detector area by us- 
ing bigger CCDs or CCD mosaics. Some dedicated IFSs 
exist but often a multi-purpose spectrograph can be given 
3D capability by inserting an integral-field unit (IFU) into 
the focal plane as if it were a standard slit mask. 

The first and simplest IFSs were derivatives of the early 
fiber multi-object spectrographs (MOS). In this design 
the input ends of the optical fibers are packed as closely 
together as possible in the focal plan and then the fibers 
reformat the 2D field into a pseudo slit at the entrance to 
the spectrograph ( [Vanderriest| [l980|. Since these IFUs 
are relatively easy to construct and incorporate with ex- 
isting spectrographs quite a few have seen operation. 



amples of productive lenslet-only IFSs are TIGRE on the 
CFHT (|Baconetal.||T995]l, OASIS and SAURON on flie 



including SILFID at CFHT (Vanderriest & Lemonnier 
[1988), INTEGRAL at the WHT ( |Arribas et al.| |I998ll 



and DensPak and Spa rsPak at WIYN (Bar den & Wade 
[T9881 [Bershady etlE] [2004) . The disadvantages of tiiis 
design are that the spatial coverage is not contiguous be- 
cause of the cladding needed to support the fibers and 
that throughput losses occur due to focal ratio degrada- 
tion (FRD) caused by feeding the fibers with a slow tele- 
scope beam (see AlUngton-Smith et al. 2002 1. 



The problems of the fiber-only IFU can be avoided by 
using a microlense array to sample the focal plane (Cour-| 



tes 1982 1. The microlenses give full, contiguous spatial 



coverage and focus an image of the telescope pupil on the 
spectrograph entrance. In lenslet-only IFSs (Figure[2]l the 
spectrograph input is an array of micropupil images. If 
grisms are used to disperse the light then the optical de- 
sign can be quite simple and have high throughput. How- 
ever, a blocking filter must be used to limit the spectral 
coverage in order to avoid spectral overlap in the dis- 
persion direction. Thus, the spectral coverage is always 
limited. Also, in the cross-dispersion direction neighbor- 
ing spectra overlap at different wavelengths. The spectra 
cannot be positioned too closely together and deconvo 
lution is mandatory to extract the spectra (IBacon et al. 
pool). Field of view is made adjustable by fore-optics 
that change the scale (arcsec/mm) at the lenslet array. Ex- 



WHT (|Bacon et al. 



2001] ), and OSIRIS on Keck ( jLarkln] 



et al. 2006[ l. Figure 3J shows a classic IFS result, the kine 



matically decoupled central disk in a velocity map from 
SAURON ( |Daviesetal.|[200T] i. 

Lenslet arrays are also used to focus pupil images onto 
the entrances of fibers. Lenslet-fiber IFUs have full spa- 
tial coverage and avoid the FRD losses of fiber-only de- 
signs. Another advantage is that the spectra can be more 
tightly packed on the detector since neighboring spectra 
on are adjacent on the sky and in wavelength (AUington- 
|Smith et al.|f2002 ). This makes more efficient use of de- 
tector area and helps produce larger fields-of-view. De- 
convolution of overlapping spectra is still helpful, but not 
mandatory. With IFUs of this type is often possible to 
trade field-of-view for spectral coverage. For example, in 
the Gemini GMOS IFUs ( AUi ngton-Smifli et al., 2002| 
and the ESO VIMOS IFU (Bonnevi lle et al.| 12003) 1 die 
full field is divided into equal sections. Fibers transfer 
the light from each section to pseudo slits at the spec- 
trograph(s) entrance(s). When the largest field-of-view 
is required blocking filters must be used to limit spectral 
coverage and avoid overlap. When higher spectral resolu- 
tion or longer wavelength coverage is desired then some 
of the pseudo slits can be masked off to increase spec- 
tral coverage at the expense of FOV. This type of IFU has 
become very popular In addition to the GMOS-N/S and 
VIMOS IFUs mentioned above, other examples include 
CIRPASS on Gemini, the A AT, and the WHT ([Parry et 
[d] 12004), PMAS on Calar Alto (|R ofli et al.|| 20051l, SPI 
RAL at the AAT (Sharp'et al'7 2()06), MPFS on the SAO 
6-m, IMACS on Mag ellan (Schmoll et al 2004), and 
FLAMES on the VLT ([Pasquini et al.||2002|l'.' ' 



The final common type of IFS is the image slicer. In this 
design the focal plane of the telescope is sampled by a 
stack of diamond-turned slicing mirrors. An optics train 
then reformats the 2D field (x, y) into a ID pseudo slit at 
the spectrograph entrance. The 2D spectra (a;, A) from the 
slices are lined up side-by-side on the detector Slicer IFU 
units can be very compact and the all-reflecting optics 
make them good for cryogenic environments. Therefore, 
they are often used for near and mid-IR IFUs. Slicer IFUs 



currently in use include 3D from MPE (|Weitzel et aL 
T996|), PIES a t Palomar ([Murphy et al.| [2000| l, UI STTt 



UKIRT (Ram say Howat et al.[[2006f 7E SI at Keck ([SheT 
msi|2006), GNIRS and NIFS at Gemini ([AlUngton-Smith 
et al.| [20061 [Hart et al.| [2003] [Miller et al.| 12006) 1, and 
SPIFFI at the VLT ( [Iserlohe et al.[|2004^ 



The native data format for each of these instruments tends 
to be unique. For fiber and lenslet IFSs the reduced 
dataset is a set of ID spectra, one for each lenslet/fiber 
spaxel, and a table of the relative positions of each spaxel 
on the sky. Often the lenslets are hexagonal so this type 
of data does not naturally produce a rectangularly grid- 
ded 3D cube. The slicer IFUs produce a set of "longslit" 
spectra that are easier to place into a regular cube. 
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Figure 2. A schematic of the three main types of integral field spectrographs: lenslet-only, lenslet+fiber (or fiber only), 
and image slicer The three techniques can produce very similar 3D data but the information density on the detector 
increases from top to bottom (lenslet only to image slicer). (Figure credit: Jeremy Allington-Smith) 
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Figure 3. The integrated light (a) and velocity field (b) of the inner region of the elliptical galaxy NGC 4365 \Davies et\ 



al. 2001 1. The galaxy has a small inner disk rotating around the minor axis but stars in the outer galaxy are rotating 



around the major axis. 
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2.3. IFU Surveys 



2.5. IFUs on JWST 



IFSs are increasingly being used for large surveys. The 
SAURON instrument on the WHT was conceived as 
a survey instrument and the initial survey of about 70 
galaxies has already revealed important new insight into 
the dynamics and formation of elliptical galaxies (Em- 
Isellem et al.| [2007l |Cappella ri et al.|[2007l l. Atlas3D is an 
extension of the SAURON survey that will observe nearly 
200 early-type elliptical galaxies to produce a complete 
sample of early-type galaxies brighter than Mk ~ —21.8 
within 40 Mpc ( |Krajnovi"3 2001) . New IFSs for surveys 
are now under development. VIRUS on the HEX is a 
massively-modular IFS being constructed for the HET- 
DEX galaxy redshift and dark energy experiment ( [Hill et| 



al. 2004 1. KMOS is a NIR IFS being built for the VLT 



which will have 24 sheer units with pickoff mirrors on 
deployable arms (Sharpies et al. ^2006). It will compile 



statistics on the kinematics of distant galaxies similar to 
what SAURON has done for nearby galaxies. Another 
new IFS for the VLT is MUSE ( jBacon etaLl[2004 ). It will 
have 24 spectrograph modules making up a 1 arcminute 
field. It is designed to work behind a ground-layer adap- 
tive optics system that should routinely produce 0.3 — 0.4 
arcsec image quality. It will be a general-purpose IFU 
studying deep fields of distant galaxies, nearby galaxy 
kinematics, and solar system objects. 



2.4. IFUs with Adaptive Optics 



If the 1990s was the era when the 8-10 meter became 
the norm, the 2000s is the time of the common use of 
adaptive optics (AO) and laser guide stars (LGS). Such 
systems are in use or currently under development at 
Keck, Gemini, ESO and the GTC they will probably be 
an integral part of the next generation of giant telescopes 
( |Crampton & Simard| |2006i Cun ningham et al.| |2006| l. 
These systems produce near-diffraction-limited images 
that rival or surpass the images from smaller space tele- 
scopes. To take advantage of the image quality spec- 
trographs with AO should have correspondingly narrow 
slits. However, normal variations in image quality and 
the "breathing", or plate scale changes normal to single- 
star AO systems could lead to significant slit losses with 
a normal slit spectrograph IFSs overcome these problems 
since slit losses are not an issue. IFSs used with AO sys- 
tems have small fields, normally around 3 arcseconds, be- 
cause of the fine sampling needed, but they provide spec- 
troscopy while preserving the high image quality. IFSs 
being using with AO systems now include OASIS at the 
WHT, GRIF at CFHT, NIFS at Gemini North, SINFONI 
at the VLT, and OSIRIS at Keck. An IFU will also be 
used with the high-contrast AO system of the Gemini 
Planet Imager (Macintosh et al. 2006j). 



IFUs will also be major components of the instrumen- 
tation for the 6.5-meter James Webb Space Telescope 
(JWST) being constructed by NASA, ESA, and the CSA. 
JWST will include three 3D spectroscopy technologies. 
A Fabry-Perot tunable filter to be provided by the CSA 
with a 2.2 arcminute field of view will be part of the 
fine guidance sensor (Rowlands et al. 2004) l. It will 



provide narrow-band imaging with R 100 and will also 
have a chonographic mode. Also, the ESA NIRSpec 
spectrograph will contain a sheer IFU with a 3 arcsec- 
ond field-of-view and 0.1 arcsec shces ( |Arribas et al.| 
|2005[ ). Finally, the medium resolution spectrometer in 
the mid-infrared instrument MIRJ is a sheer IFS with four 
slicers in separate wavelength channels split by dichroics 
( [Wright etall[2004l l. 



3D SPECTROSCOPY AND THE VIRTUAL OB- 
SERVATORY 



The increasing popularity of optical/IR 3D data and the 
large volume of 3D data that will soon be in archives 
makes it very important that 3D spectroscopy is sup- 
ported by the VO. As with any form of astronomical data, 
the main issues that have to be addressed are data descrip- 
tion, storage format, discovery, data transmission, and 
analysis tools. 3D data provides an interesting challenge 
because it is a mixture of imaging and spectroscopy. For- 
tunately, this means that existing protocols and tools for 
handling imaging and spectroscopic data separately can 
serve as a basis for handling 3D data. A goal is that an 
astronomer could be able to use the VO to find and com- 
bine multi-wavelength 3D data on an object (a 3D SED) 
and then do analysis on this meta-cube. With the addition 
of multiple observations to give the time information, we 
are approaching the goal of the fully general observation. 



3.1. Data Description 

A meta-data structure, or general data model, for de- 
scribing optical/IR 3D data has been developed using the 
IVOA Simple Spectra Access Protocol (Chilingarian et 
[aL]|2006| ). More details about this standard and an exam 
pie apphcation are given in (Chilingarian et al. 2007 1. 



3.2. Storage Format 

The storage format for 3D data is an important consid- 
eration for discovery, transmission, and analysis. As has 
been mentioned, each IFU tends to have its own data for- 
mat. This is partially due to the fact that the most com- 
monly used data reduction packages such as IRAF, MI- 
DAS, and IDL, do not have a common data structure for 
3D data. Also, there is the problem that different types of 
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instruments produce different types of data: stacks of im- 
ages; one spectrum per spaxel; or a series of 2D "longslit" 
spectra. A common data format for reduced data would 
help promote the use IFU use by making it easier to de- 
velop analysis tools. A common format would also sim- 
plify making IFU data VO compliant. 

One option for a standard IFU data format is the Euro3D 
format developed by the European OPTICON Euro3D re- 
search training network (Kissler-Patig et al.) |2004| . The 
format is a 3D FITS table in which each row contains the 
information from one spaxel. Columns can contain co- 
ordinate information and data, data quality, and statistics 
vectors. It is general to all IFS instruments after the in- 
strument signature is removed. It also has a very general 
spatial description of a spaxel that does not require re- 
sampling to a 3D linear cube and that supports adaptive 
binning ( [Cappellari & Copin 2003) 1. 



Another option for a common format is the linearly sam- 
pled 3D FITS cube. This gives the data a similar format 
to radio and 3D spectrometer data cubes and it makes the 
data a bit easier to visualize with common display tools. 
However, it would require that many IFS datasets be re- 
sampled (which can often degrade data) and it doesn't 
lend itself to adaptive binning. Data that is not on a 
regular grid should not be forced into one unnecessarily. 
However, because of the meta-data description it should 
be possible to support both Euro3D and datacube for- 
mats. 



3.5. Visualization and Analysis Tools 



As has been mentioned several times, visualization and 
analysis of 3D data is challenging because of the com- 
bination of needs for manipulating images and spectra. 
It is always important to be able to treat a 3D dataset as 
a set of images. It is often useful to loop through the 
planes like a movie. For IFS data this means that the 
tool needs to reconstruct the images and show any bin- 
ning that was done. Registering and either overlaying or 
overplotting different planes or different images are also 
needed. Also, it should be possible to select a wavelength 
range and inspect the resulting 2D spatial map. Inversely, 
one should be able to select a region from a map and then 
inspect either the individual spectra or the summed or av- 
eraged spectrum. As with imaging, overplotting is an im- 
portant feature. Two 3D visualization tools that do many 
of these things are the Euro3D tool (Sanchez, 2004) and 
QFitsView (Abuter et al., 2006). However, neither is VO 
compliant at this time. Another potentially powerful vi- 
sualization option is to treat the spatial and spectral di- 
mensions together A few visualization labs and authors 
have been experimenting with generating 3D views of 
IFU data by using volume rendering, stereoscopic views 
or "3D" display monitors, and holographic imaging (see 
Figure 4|l. These techniques make full use of the data and 
can help astronomers draw the correct physical interpre- 
tation. Generating good visualization tools for 3D data is 
a challenge but also an important opportunity for the VO 
community. 



3.3. Discovery 

Data discovery needs to answer the user questions "What 
data is available?" and "Is it useful to me?". If the answer 
to the second question is "yes" then the astronomer will 
continue to work on the data by downloading it or analyz- 
ing it with VO tools. For 3D data good query interfaces 
must be joined with both imaging and spectral visualiza- 
tion tools. Searches need to be done both in terms of 
position and wavelength. It should be possible to visual- 
ize the footprint of IFU data on a digital sky image. Then 
it should be possible to preview both image and spectral 
representations of the 3D dataset. A first-generation tool 
that fulfills many of these requirements is described by 
ChiUngarian et al. ( |2007| l in this volume. 



Until now we have been treating the spectral and spatial 
domains independently, however, certain operations on 
3D data require that they be treated simultaneously. One 
example is the need to subtract a spectral continuum from 
a region using an annulus around it. This would involve 
fitting in both wavelength and position. Also, both spatial 
and spectral PSFs change with wavelength so a deconvo- 
lution of either quantity involves deconvolution in both. 
For example, atmospheric refraction is a shift of spatial 
position with wavelength and correcting for it requires 
operations in all dimensions of the dataset (Arribas et al.] 
|I99 9). Finally, adaptive binning is a means of combining 
adjacent spaxels until the spectral signal-to-noise is suffi- 
cient for a particular analysis. All these operations mean 
that the errors in the spectra of adjacent spaxels become 
correlated so error propagation is very important for the 
final assessment of signal-to-noise. 



3.4. Data Transmission 

Full Euro3D or datacube files can be many megabytes 
in size and new instruments like VIRUS and MUSE will 
produce vastly larger datasets than current instruments. 
Therefore, it is important that it be possible to extract re- 
gions of a dataset. Because of the size and complexity of 
3D data, data providers should experiment with provid- 
ing processing and analysis on the server. However, the 
users should always have the option to download the data 
for local analysis. 



The analysis needs in only the spectral dimension are 
fairly standard, only it must be easy to apply the same 
analysis to all spaxels in a dataset. The main needs are 
single and multiple line fitting, radial velocity and veloc- 
ity dispersion measurements using observed templates or 
models, and the measurement of line fluxes, equivalent 
widths, and line indices. Radial velocity measurements 
and stellar template fitting to separate the emission line 
and absorption line spectra are good applications for the 
VO since the template observations or models could be 
accessed via the VO. This also meshes with the VO focus 
on SED fitting which should be generalized to 3D. 
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Figure 4. Volume renderings of 3D data cubes, a) Lya 
emission and absorption from a z = 3.09 galaxy ob- 
served by SAURON ( Wilman et al. 2005 1, b) Represen- 
tations of the H2 and Br"f emission from the planetary 
nebula NGC 7027 ( Cox et aL\ 2002). The outflow cavi- 
ties sketched in Figure\ljare clearly seen. 



4. SUMMARY 



Three dimensional spectroscopy techniques in the optical 
and NIR are a relatively new but rapidly growing in im- 
portance. With 3D instruments on most large telescopes, 
several large 3D survey instruments under development, 
and 3D instruments being prepared for JWST, there will 
be en ever growing amount of 3D data in archives which 
future investigators will need to access easily in order to 
take full advantage of them. Therefore, the VO needs 
to support the discovery, transmission, and visualization 
and analysis of 3D (and higher dimensional) data. There 
are challenges to handling these more general datasets 
because of the union of spatial, spectral, and eventually 
time series and polarization analysis. However, this is 
also an important opportunity because tools for handling 
3D data are not well developed. Good tools would draw 
people with IFU data to the VO and it would allow the 
joining of 3D datasets from all wavelength ranges. The 
discovery and exploration would allow for the creation of 
the most complete combined datasets possible and of data 
products such as the spatially resolved SED. The capabil- 
ity to handle such datasets will encourage data providers 
to make their data VO accessible and give investigators 
more incentive to use the services, thereby maximizing 
the potential for scientific discovery. 
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Once a 2D map is produced it should be possible to 
perform all standard imaging analysis on it. This in- 
cludes calculating statistics, aperture photometry, ellipse 
or isophote fitting, fitting surfaces, and the subtraction of 
background regions. In addition, many maps may not 
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